Geometric reheating after inflation 
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Inflationary reheating via resonant production of non-minimally coupled scalar particles with only 
gravitational coupling is shown to be extremely strong, exhibiting a negative coupling instability 
for £ < and a wide resonance decay for ( > 1, Since non-minimal fields are generic after 
renormalization in curved spacetime, this offers a new paradigm in reheating - one which naturally 
allows for efficient production of the massive bosons needed for GUT baryogenesis. We also show 
that both vector and tensor fields are produced resonantly during reheating, extending the previously 
known correspondences between bosonic fields of different spin during preheating. 

PACS: 98.80.Cq, 04.62,+v 



3 



> 

o> 
o 

On 



i 

Oh! 
CD . 
43 ■ 



X 



I. INTRODUCTION 

Typical modern incarnations of inflation arise within 
supergravity, string or GUT theories where the inflaton, 
4>, is only one of many fields. Studies of inflation in- 
cluding couplings to these other fields, as required to re- 
heat the universe after inflation, yield extremely complex 
dynamics M and are little investigated beyond hybrid 
models 0. Here we turn to a minimalist view in which 
preheating occurs with the inflaton coupled only 

gravitationaly to other fields. We call this geometric re- 
heating to emphasize its gravitational origin. The most 
powerful example of this new mechanis m is provided by 
non-minimally coupled fields (section II B ) , where the 
strength of the effect is due to the congruence of two 
facts: (1) the Ricci curvature oscillates during preheating 
and (2) the non-minimal coupling, £, is a free parameter. 
The first ensures that there is resonance, the second that 
the effect is non-perturbative. 

The possible importance of this effect is motivated 
by renormalization group studies in curved-spacetime 
|^| [Tof , which have shown that even if the bare coupling, 
£0 is minimal, after renormalization £ =/= gcnerically. 
Since we are particularly interested in the preheating 
realm which occurs when inflation ends near the Planck 
scale, we are near the ultra-violet (UV) fixed points of the 
renormalization group equations. While the UV fixed 
points may correspond to a conformally invariant field 
(to = 0, £ = jjj), in different GUT models the coupling 
may also diverge, |£| — ► oo, in the UV limit (||ll]]. ^ n 
both cases the nature of geometric reheating is very dif- 
ferent from the standard models based on explicit self- 
interactions or particle-physics couplings between fields 
(see e.g. |). ' 

Further we shall study the gravitational production 
of spin 0,1 and 2 particles due to the expansion of the 
universe during preheating, and will show that a unified 
treatment in terms of parametric resonance exists. This 
is shown by reducing the evolution equations to general- 
ized Mathieu form: 



[A(k) - 2qcos2z]x = 0, 



(1) 



with time-dependent parameters [^. 

The Mathieu equation has rapidly growing solutions 
controlled by the Floquet index In the case that 
1 3> q > 0, the Floquet index in the first resonance 
band of the Mathieu equation is given by /ife = ((q/2) 2 — 
(2k /m - l) 2 ) 1 / 2 IUI. This can be extended to give /<tf 
in the N-th resonance band as long as A > and 
2iV 3 / 2 > q: 
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where 5 varies in the interval [— 7t/2,0] and fi^ -C 1. 

When A(k) < a qualitative change occurs and the 
dominant effect comes from the fact that one effectively 
has an inverted harmonic oscillator yielding the negative 
coupling instability jl4j|. In this case the Floquet index 
can be as large as fik ~ kl 1 ^ 2 , there are no stability 
bands to speak of and the typical variances are larger by 
a factor \q\ x / 2 than in the A(k) > case. 

To be concrete, consider the case of a scalar field 
in a FLRW universe (g M „ = diag(— 1, a 2 (t)/(l — 
Kr 2 ),a 2 (t)r 2 ,a 2 (t)r 2 sin 2 6) K = ±1,0 is the curvature 
constant). We shall restrict^] ourselves to the quadratic 
potential, 



* While it is known that the Mathieu formulation is insuffi- 
cient H in some respects, and has lead to the introduction 
of other approximations - principally that of stochastic reso- 
nance |3| - the Mathieu equation remains a powerful diagnos- 
tic test for the strength of particle production. 

'We note that application of stochastic resonance methods 
to the vector, tensor and non-minimal scalar fields of this 
paper for the potential V = j4> 4 requires an extension of the 
existing theory to scattering in a quartic potential as opposed 
to the standard quadratic potential ^Q. 
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For K = 0, the latter potential gives an oscillatory be- 
haviour of the field, <fi = $ sm(m#t), with <f> ~ l/m</,t. 
In what follows we shall try to preserve maximal gener- 
ality; we denote with = results which are derived specif- 
ically for the potential . We use natural units with 
k = 8tt,G= 1. 

The energy density and pressure for a minimally cou- 
pled scalar field, treated as a perfect fluid, are \i = 
+ V((f>)), p = - V{<f>)). This breaks down 

if the field is non-minimally coupled (an imperfect fluid 
treatment must be used), if the effective potential is not 
adequate Q, or if large density gradients exist. The 
FLRW Ricci tensor is CI 
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where i, j — 1..3. The Ricci scalar is: 
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The Raychaudhuri equation for the evolution of the ex- 
pansion O = 3d/a is [j given by: 



3K 

^2" 



while the Friedmann equation is 



6 2 + -5- 



3 KA t = 3k ( + y(</>) 



(6) 



(7) 



As an example, when K = and d/(am ( p) <C 1, one 
may solve Eq. (j?]) perturbatively 0, to obtain: 
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to first order in a/{am$). This is only valid after preheat- 
ing when $ <C 1 but shows that the expansion oscillates 
about the mean Einstein-de Sitter (eds) pressure-free so- 
lution. Eq. M) can be integrated to give the scale factor: 



a(t) = a exp 



nn2m c j > t 2ci(2?7j0i) 



(9) 



where a — i 2 / 3 is the background EDS evolution, and 
ci(m^,t) — — J cos(m<pz)/ zdz. This example explicitly 



'•The expansion is generally defined as O = u a - a where u a is 
the 4- velocity and ; denotes covariant derivative fll6| |. 



demonstrates how temporal averaging (which yields a) 
removes the resonance. 

Via Eq.'s (^,0) one can systematically replace all fac- 
tors of d, d with factors of (/> and V(<p) terms^]. In this way 
one can show that the vector and tensor wave equations 
take the form of Mathieu equations during reheating Q . 



II. SCALAR FIELDS 

Consider now the effective potential: 

N N 

2 m ^ + 2 



1 N 1 N 

Xu) = v(4>) + mlxl + 7, E > (io) 



describing the inflaton with potential V(<f>) coupled only 
via gravity to N scalar fields, which have no self- 
interactions, masses m„ and non- minimal couplings £„. 
The equation of motion for modes of the v-th. field is: 



r k + QXk+[y+rni+^Rjxk = 0, (11) 

From Eq (||,|||j]) the Ricci scalar is given by |**| : 

R = -k^? + 4nV{(j)) (K = 0). (12) 

A. The minimally coupled case 

Consider = 0. Then ( pd| ) reduces, on using Eq's 
(1,0), to: 



d 2 (a 3 / 2 X k) , / k 



dt 2 
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There exists parametric resonance because the expansion 
O oscillates. The potential (||) yields equation (|l|) (K — 
0) with time-dependent parameters: 



A(k,t) 



k 
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(14) 



From this we see that the production of particles is re- 
duced as m„ increases. Indeed, since A — ► m^/mi, q — > 
due to the expansion, production of minimally coupled 



^Indeed, a useful combination is 2Q + G 



!k7710$ 2 cos(2m^,t) 



-3np 



**Assuming that at the start of reheating the inflaton is the 
dominant contributor to the energy density of the universe. 
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bosons is rather weak and shuts off quickly due to hori- 
zontal motion on the stability chart. We stress that the 
production is, however, much stronger than that obtained 
in previous studies where the scalar factor evolves mono- 
tonically |]l9| . This mild situation changes dramatically 
when a non-minimal coupling is introduced. 



B. Non-minimal preheating 

Now include the arbitrary non- minimal coupling £„. 
Using Eq. jl2] ) one can reduce Eq. (11) to (K — 0): 



dt 2 
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« (« 3/2 Xfe)=0. (15) 



Defining a new variable z — m,pt+7r/2, Eq. (15) takes the 
form of equation ([!]) with time-dependent parameters: 



A(k,t) 
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The crucial observation is that since £ is initially free 
to take on any value A(k) is neither restricted to be 
positive nor small. 

From Eq. ( [l6| ) it is clear that A{k) < for sufficiently 
negative £. The posibility of negative A was the thesis of 
the work by Greene et al |1J]. However, in their model, 
this powerful negative coupling instability was only par- 
tially effective due to the non-zero vacuum expectation 
value acquired by the x field due to its coupling, g, with 
the inflaton. Here we only have gravitational couplings 
and the same constraint is removed. 

Negative A (induced when q < 0) implies that the 
physical region of the (A, q) plane is altered. Instead 
of A > 2\q\ we have A > tt$ 2 - 2|g|/3. Now when 
2\q\/3 > \A\ > 1 we have ^ k ~ \q\^ 2 ~ (6Tr\^\)^ 2 ^> along 
the physical separatrix A = ir$ 2 — 2|g|/3. Since the renor- 
malized |£| may have very large values, this opens the 
way to exceptionally efficient reheating - see Figs. fflJa) - 
via resonant production of highly non-minimally coupled 
fields with important consequences for GUT baryogenesis 
p3 and non-thermal symmetry restoration. 



^The only constraints that one might impose are that the 
effective potential should be bounded from below and that the 
strong energy condition, R a bU a u b > <S> T a bU a u b > — T/2, 
be satisfied. The first is difficult to impose since R oscillates 
and the second since one should use the renormalized stress- 
tensor, (Tab)- 




FIG. 1. The evolution of the k — mode (m^/m| ~ 1), 
as a function of time and the non-minimal coupling param- 
eter £. For positive £ the evolution is qualitatively that of 
the standard preheating with resonance bands. However, for 
negative A (negative f) the solution changes qualitatively and 
there is a negative coupling instability. There are generically 
no stable bands and the Floquet index corresponding to — 1£| 
is much larger, scaling as ~ |£| 1//2 - 




FIG. 2. A slice of the spectrum in fig. (1) at t — 5 as 
a function of the non-minimal coupling £. The qualitative 
differences between £ < and £ > are clear. 

For example, let us consider ~ 2x 10 13 GeU as 
required to match CMB anisotropies AT/T ~ 10~ 5 . 
Then GUT baryogenesis with massive bosons x with 
m x > 10 14 GeU simply requires £ < — (7r$ 2 ) -1 , with $ 
in units of the Planck energy. Instead if one requires 
the production of GUT-scale gauge bosons with masses 
rrigi, ~ 10 16 GeV this is still possible if the associated 
non-minimal coupling is of order £ ~ — 10 3 . Such cou- 
pling values have been considered in e.g. p0[ . The mas- 
sive bosons with m x ~ 10 14 GeV can be produced in the 
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usual manner via parametric resonance if £ > 0, but this 
process is weaker (c.f. [fjl[). 

Since the coupling between <f> and x is purely gravi- 
tational, backreaction effects in the standard sense (see 
H[|) cannot shut off the resonance. The inflaton con- 
tinues to oscillate and produce non-minimally coupled 
particles, receiving no corrections to m\ e ^ from (<5x 2 )- 

To estimate the maximum variance (5x 2 ) is therefore 
rather difficult. The standard method is to establish the 
time when the resonance is shut off by the growth of A(k) 
which pushes the k — mode out of the dominant first 
resonance band. For this we must understand how A(k) 
changes as the %-field gains energy and alters the Ricci 
curvature. If we assume that most of the energy goes 
into the xo mode, justified in the £ < case then the 
change to the Ricci curvature is 5R X = 8w(E — S), where 



E = G, 



If 



m 2 x xl 



12£xoXo- 
a 



(17) 



is the T uu component of the \ stress tensor, 



8 = 



3G 



eff 



1 + l92TrG eff ex 2 o 



±1 rn\xl 



*2„,2 ; 



< ( ^ - XoXo^ - rn 2 x xl ) + M»C\r,E 



(18) 



is the spatial trace of the stress tensor T% corresponding 
to 3p in the perfect fluid case and G e // = (l + Wn^x 2 )" 1 
is the effective gravitational constant. Since xo is rapidly 
growing, the major contribution of SR X will be to A(k), 
causing a rapid vertical movement on the instability 
chart. Once SA + A > 2\q\ + |(j , | 1/ ' 2 , the resonance is shut- 
off. If £ < 0, most of the decaying <fi energy is pumped 
into the small k modes (see fig. |2|). Subsequently we 
expect the oscillations in xo to produce a secondary res- 
onance due to the self-interaction and non-linearity of 
Eq's (0,©. 

The case of a A x x 4 /4 self-interaction provides another 
mechanism that may be dominant in ending the reso- 



and hence A(k), gains 
■ (Sx 2 ) which shuts off the 



nance: namely e j ^ , 
tions proportional to A Al 

onance jlij leaving a peak variance of order (Sx 2 ) 



res- 



')/A x (assuming that £ < 0). If £ > the 



variance is smaller by a factor \q 



1/2 



III. THE VECTOR CASE 

Until now, reheating studies have been limited to 
minimally-coupled scalar fields, fermions and gauge 



bosons p3[ . In the case of vector fields the minimum 
one can do to preserve gauge-invariance is to couple to 
a complex scalar field via the current since real scalar 
fields carry no quantum numbers. We consider here only 
vacuum vector resonances, however. 

A massive spin-1 vector field in curved spacetime sat- 
isfies the equations: 



(-v a v Q 



0A b + R\A a - 



(19) 



These equations are equivalent to the Maxwell-Proca 
equations for the vector potential A a only after an appro- 
priate gauge choice which removes one unphysical polar- 
ization. In our case we shall use the so-called tridimen- 
sional transversal constraint: 



An = 0, V*A = 



(20) 



This set is equivalent to the Lorentz gauge, although 
it doesn't conserve the covariant form of the latter. 
Nonetheless, in either case, gauge-invariant quantities 
such as the radiation energy density, are unaffected. 

In a FLRW background, the Ricci tensor is diagonal, 
which together with the gauge choice (20) and expansion 
over eigenf unctions, ensures the decoupling of the set of 
equations (0). We can reduce the system to a set of 
decoupled Mathieu equations. The Ricci tensor is (see 
Eq. |): 



(21) 



which leads to the Mathieu parameters for for the spatial 
components (a 3 ^ 2 A l ): 



A(k) 



9 2 



-4 + 2q 



K$ 2 , , 

-5" ( 22 ) 



showing that vector fields are also parametrically ampli- 
fied, albeit weakly, during reheating as in the scalar case. 



IV. THE GRAVITON CASE 

It has been shown using the electric and magnetic parts 
of the Weyl tensor |lj| that there exists a formal anal- 
ogy between the scalar field and graviton cases during 
resonant reheating. Here we will show that the corre- 
spondence also holds in the Bardeen formalism. The 
gauge-invariant (at first order) transverse-traceless (TT) 
metric perturbations /ly describe gravitational waves in 
the classical limit. In the Heisenberg picture one expands 
over eigenfunctions, Y a b of the tensor Laplace-Beltrami 
operator with scalar mode functions hu, which satisfy the 
equation of motion: 



hk + &h k + 



2K 



h k = 0, 



(23) 



"In the case £ 3> 1 one needs to use (5x 2 ) instead. 



or equivalcntly 
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(a 3/2 h k y+ + b) (« 3/2 ^) = , (24) 

where p = K,(<j) 2 /2 — V) is the pressure. This gives a 
time-dependent Mathieu equation (c.f. Eq. [l4]) with pa- 
rameters: 



A(k) 



q = 



16 



(25) 



In this case, a negative coupling instability is impossi- 
ble and only for <J> ~ M p i is there significant graviton 
production. Note, however, that if temporal averaging is 
used, the average equation of state is that of dust, p = 0. 
Eq. (|24|) then predicts (falsely) that there is no reso- 
nant amplification of gravitational waves since the value 
of q corresponding to the temporarily averaged evolution 
vanishes. 



V. CONCLUSIONS 

We have described a new - geometric - reheating chan- 
nel after inflation, one which occurs solely due to grav- 
itational couplings. While this is not very strong in 
the gravitational wave and minimally coupled scalar field 
cases, it can be very powerful in the non-minimally cou- 
pled case, either due to broad-resonance (£ 3> 1) or nega- 
tive coupling (£ < 0) instabilities. Particularly in the lat- 
ter case, it is possible to produce large numbers of bosons 
which arc significantly more massive than the inflaton, as 
required for GUT baryogenesis. It further gives rise to 
the possibility that the post-inflationary universe may 
be dominated by non-minimally coupled fields. These 
must be treated as imperfect fluids which would thus al- 
ter both density perturbation and background spacetime 
evolution, which are known to be significantly different 
pij than in the simple perfect fluid case. We have fur- 
ther presented a unified approach to resonant production 
of vector and tensor fields during reheating in analogy to 
the scalar case. 

Future work should examine in greater detail the pos- 
sibility of negative coupling instability in the vector field, 
backreaction issues in the non-minimal case, and the sit- 
uation in potentials with self-interaction relevant to sym- 
metry restoration. 
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